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At low energy scales charge transport in the insulating 
Si:P is dominated by activated hopping between the lo- 
calized donor electron states. Thus, theoretical models 
for a disordered system with electron-electron interac- 
tion are appropriate to interpret the electric conductivity 
spectra. With a newly developed technique we have mea- 
sured the complex broadband microwave conductivity of 
Si:P from 100 MHz to 5 GHz in a broad range of phos- 
phorus concentration n/ric from 0.56 to 0.95 relative to 
the critical value ric = 3.5 x 10^^ cm^"^ corresponding 
to the metal-insulator transition driven by doping. 



At our base temperature of T = 1 . 1 K the samples are in 
the zero-phonon regime where they show a super-linear 
frequency dependence of the conductivity indicating the 
influence of the Coulomb gap in the density of the impu- 
rity states. At higher doping n ric, an abrupt drop in 
the conductivity power law (Ji{lo) ^ uj°' is observed. 
The dielectric function ei increases upon doping fol- 
lowing a power law in (1 — n/ric)- Dynamic response 
at elevated temperatures has also been investigated. 
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1 Introduction In a material of such an industrial im- 
portance like Si:P there are still questions open about the 
low-energy excitations from the ground state. In particu- 
lar, the influence of electron-electron interactions on the 
hopping transport and the critical behavior at the metal- 
insulator transition (MIT) have attracted much attention 
since decades. Though theoretical predictions on the T-^O 
frequency-dependent response of an interacting disordered 
system have already existed since many years H], experi- 
mental data still remain scarce and results obtained by dif- 
ferent groups and in different parameter ranges lack con- 
sistency. Especially in the microwave range, from tens of 
MHz till tens of GHz, with wavelengths from a few mil- 
limeters to a few meters, no better means than the reso- 
nator techniques have been used for a long time to study 
the broadband dynamic conductivity of doped semicon- 
ductors. With a novel and advanced method of measuring 
the broadband (100 MHz to 5 GHz) complex microwave 
conductivity of that material class, we focus on the con- 
ductivity power law in Si:P at T^O and at elevated tem- 
peratures in a broad range of donor concentration close to 



the MIT. 

At concentrations of phosphorus in silicon below the 
critical value of Uc = 3.5 x 10^^ cm^'^, the donor electron 
states are strongly localized due to disorder in Anderson 
sense |2l and the corresponding wave functions resemble 
those of a hydrogen atom t3 |. Since some degree of com- 
pensation by impurities of the opposite type is considered 
inevitable, charge transport at low excitation energies is by 
variable-range hopping between the donor sites, randomly 
distributed in space IS]!?]. 

1.1 Dynamic conductivity The static conductivity 
(Tdc of the insulating Si:P vanishes when T ^ 0. The main 
issue we address is that of power laws of the frequency- 
dependent conductivity in the variable-range hopping regi- 
me at zero temperature: 



fTl (w) 



a = ai + ia2 ■ 



(1) 



The theory of resonant photon absorption by pairs of states, 
one of which is occupied by an electron and the other one 
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is empty, yields distinct limiting results for the conducti- 
vity power law in cases, where one of the relevant energy 
scales of the problem dominates over the others: 

For the non-interacting system, where the photon ener- 
gy prevails over the electronic correlations, Mott has pre- 
dicted a sub-quadratic frequency-dependence of the con- 
ductivity ID: 

ai{uj) = -qe^ ar'^hhj'^ , (2) 

= aln(2/o/?iijj) . 

Here, 77 is a numerical coefficient lUlU, g is the state den- 
sity at the Fermi level, a is the localization radius, is the 
most probable hopping distance and Iq is the pre-exponen- 
tial factor of the overlap integral between the localized elec- 
tron states. 

Taking into account the Coulomb repulsion U{r^) if 
both states in a pair would be occupied by an electron, 
Shklovskii and Efros corrected the picture of the energy 
levels before and after a photon is absorbed and derived 
ai{uj) to be a sub-linear function of frequency, as long as 
the Coulomb interaction term dominates over the photon 
energy HI: 

ai{uj) = rie'^g'^ar^u![huj + U{r^)], (3) 



At higher frequencies, in the opposite limit, the sub-quadratic 
behavior known from Mott ^ is recovered. 

The formula ^ is derived under the assumption of an ef- 
fectively constant density of states g near the Fermi level. 
This implies another characteristic energy of the system, 
the Coulomb gap A, to be small compared to U{rS) and 
thus of no significant effect on the states participating in the 
hopping transport. The parabolic Coulomb gap forms in 
the density of the impurity states around the Fermi level as 
another consequence of the electronic correlations 13]. For 
the conductivity of interacting electrons where the Coulomb 
term U {r^) dominates over the photon energy but falls in- 
side the Coulomb gap, the reduction of the density of states 
leads to a stronger, slightly super- linear power law HJ: 

cri(w) ~w/ln(2/o/?kj). (4) 

1.2 Dielectric function It is an additional advantage 
of a phase sensitive measurement to gain the dielectric func- 
tion £1 from the imaginary part of the complex conductivi- 
ty ct = cti + i(72 H): 

a 

e = £i+«e2 = l + i . (5) 

We denote by e the full complex dielectric function of Si:P, 
relative to the permittivity Eq of vacuum, and use the SI 
units throughout this text. As the MIT is approached upon 
doping n, the localization radius a diverges according to a 
power law ||3] : 

a ~ |1 — n/nc\~^ ■ 



As a consequence, the electronic contribution to the dielec- 
tric function is also expected to diverge following a power 
law when the MIT is approached 16J : 

£1 - £si |1 - n/n-cT'' , (6) 

where £si = 11.7 is the dielectric constant of the host ma- 
terial Si. 

2 Experiments 

2.1 Samples Si:P-samples for this study were cut 
from a Czochralski-grown cylindrical crystal IT), nomi- 
nally uncompensated, with a phosphorus concentration gra- 
dient along its axis. To remove distorted surface layers [8] 
121, the samples were chemically and mechanically treated 
by well established procedures. The donor concentration 
was determined from the room-temperature resistivity ifTOl 
employing a commercial four-probe measurement system 
(FPP 5000 by Veeco Instruments). For high doping levels 
the resistivity ratio p(4.2K)//9(300K) (determined from stan- 
dard dc measurements) is consistent with Ref. ifTTI . In the 
present work, measurements on Si:P-samples with relative 
donor concentration n/ric ranging from 0.56 to 0.90 rela- 
tive to the critical value at the MIT are discussed, as sum- 
marized in Tab.[T] 

Table 1 Dopant concentration n and the value of the 
room-temperature dc resistivity pdc of the presented Si:P- 
samples. The second column lists the relative concentra- 
tion with respect to the critical value ric = 3.5 x 10^^ 
cm~'^ at the MIT. The power a of the conductivity power 
law di (uj) uj°' and the full relative dielectric constant £1 
are also displayed in Fig|4]and Fig|7] 
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n/ric 


Pdc(300K) 
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El 


(10^* cm^^) 




(SI cm) 






1.97 


0.56 


0.0162 


1.13 


23 


2.29 


0.65 


0.0149 


1.16 


24.5 


2.57 


0.73 


0.0139 


1.10 


28 


2.91 


0.83 


0.0130 


1.08 


34 


3.04 


0.87 


0.0127 


1.05 


41 


3.14 


0.90 


0.0124 


1.04 


46 
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Figure 1 Scematic of a broadband microwave spectrometer that 
employs a vector network analyzer. A solid sample terminates the 
otherwise open-ended coaxial line (after 1131 ). 



2.2 Broadband microwave spectroscopy The 

broadband microwave spectrometer, previously employed 
for measuring metallic samples 11311141 . has been upgraded 
lfT2l with microwave components, calibration technique and 
evaluation procedure suitable to study the hopping trans- 
port in doped semiconductors at temperature as low as pos- 
sible. 

The sample with a flat surface terminates the 1 .2m-long 
coaxial line, Figd] needed to transport the microwave sig- 
nal from the source to the bottom of a "'He-pumped cryo- 
stat. The reflected signal, containing all the useful informa- 
tion about the sample, returns the same way back to the test 
set of a HP 8510 vector network analyzer. The measure- 
ment is phase sensitive and yields the complex reflection 
coefficient. After an extensive calibration procedure requi- 
ring three independent low temperature measurements of 
calibration samples under reproduced conditions, the re- 
flection coefficient 5ii at the sample surface is obtained. 
The complex impedance Z is directly related to it via lil5il : 

Z-Zo 



Su = 



(7) 



Z + Z„ 

with the characteristic impedance Zq of the microwave line. 

In a next step, the complex conductivity a has to be 
extracted from the impedance Z. If material is insulating, 
there exists no direct solution, because the simple concepts 
known from metallic samples lfT4l fail here. In the case 
of an insulator, the surface impedance approach 0, useful 
for thick metallic samples, yields a wrong (i.e. too strong) 
frequency dependence of the conductivity cri . If the semi- 
conducting sample is treated as a thin highly-conducting 
film, the conductivity values are suspiciously large lfT6l . 

2.3 Advanced data analysis For these reasons, we 
have developed a general and rigorous solution for the prob- 
lem. The electromagnetic wave penetrates deep into a thick 
insulating sample, FiglH and forms a three-dimensional 
distribution quite different from that of a plane wave. Thus, 
we work with a solution for the field distribution in an in- 
tegral equation formulation, combined with the variational 
principle IITtI. Ascertaining, that the electromagnetic field 
strength decreases below 1 % at the depth corresponding to 
our sample thickness of 2 mm in the frequency range eva- 
luated, we forbear from taking into account the secondary 
reflections at the back side of the sample as in the full- 



Al contacts matching the probe 
2b- 




Figure 2 Sketch of the Si:P-sample with the aluminium 
contact-layer on top. 2a = 0.6 mm, 2b = 1.75 mm, correspond- 
ing to a 2.4 mm microwave adapter, c = 5 mm, d = 2 mm 



wave analysis at higher frequencies considered by Brom 
and collaborators ifTSll . The resulting formula connects the 
complex sample impedance Z to the wave vector k, the lat- 
ter containing the complex dielectric function e we search 
for (Note, that we hold to the convention cxp(i(kx — uit)) 
for the Fourier-transformed of the electromagnetic field in 
accordance with ||5]fT7l. The opposite sign convention is 
used in Ref. 1181 191 and by the network analyzer.) 



Zo 
Z 



b b TT 



nkc ln(6/a) 



cos tp ■ 



cxp(iA:r) 



dip dp dp', 



a 



r= (p2 + p'2_2pp'cos^)i/2_ (8) 

To solve this integral equation for k we make use of 
the practical ideas from |fT9ll . In a series expansion of the 
exponential function the terms higher than the third term 
are negligible below 5 GHz, integrating the rest one obtains 
a quadratic equation in e = ei + i£2 lfT2l : 



1 

Z 



—i2ujeo£ 
[ln(6/a)]2 



(9) 



with the geometrical integrals Ii = 0.9084 mm and I3 = 
—0.2100 mm'' numerically evaluated for the dimensions a 
and b of our coaxial probe as shown in Fig|2] 

2.4 Low-temperature ac measurements The fre- 
quency range from 100 MHz to 5 GHz for the reflection 
coefficient measurements according to section [Z2l contains 
300 frequency points: 200 equidistant points from 0. 1 GHz 
to 1 GHz and 100 from 1 GHz to 5 GHz. The spectra are 
taken at the base temperature I.l K of our ^He-pumped 
cryostat as well as at elevated temperatures using a tem- 
perature controller Lake Shore 340 in an automatized pro- 
cedure while heating the setup up to room temperature. All 
the reported samples prove to be in the zero-phonon regime 
at T = 1 . 1 K because they show a saturation of cti (T) and 
El (T) in the whole frequency range as T ^ 1 K. 

3 Results and discussion For each doping concen- 
tration n four to ten independent measurements have been 
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Figure 3 Typical spectra of the measured real part of 
the conductivity for Si:P with relative donor concentrations 
n/uc of 0.65, 0.83 and 0.90. The dashed lines are the fits 
by a two-parameter power law function. 




0.6 0.7 0.8 0.9 

phosphorus concentration n/n^ 

Figure 4 Mean value of the conductivity power a obtained 
from the tri ~ uj" fit. The dependence on the relative donor 
concentration n/uc becomes stronger as rt ^ = 3.5 • 



conducted. From the best spectra (lowest noise and smal- 
lest influence of standing waves) mean values of the con- 
ductivity power a and the dielectric constant ei have been 
determined. 

3.1 Frequency dependence of the conductivity in 
the zero-phonon regime In Figl3]the measured real part 
of the frequency-dependent conductivity is plotted on a 
log-log scale to identify the power law. The fits by a two- 
parameter function ai (oj) ~ const • oj" crossing the origin 
are shown by the dashed lines. In Fig|4]the mean values 
of the power a are plotted against the relative dopant con- 
centration. The frequency dependence of the conductivity 
clearly follows a super-linear power law in the whole do- 
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Figure 5 Temperature dependence of the conducti- 
vity at fixed frequency on example of relative doping 
n/ric =0.83, typical for all the investigated Si:P-samples. 



ping range, where the exponent decreases slightly with do- 
ping; this effect becomes even stronger when n is increased 
further. From this we infer, that hopping transport takes 
place deep inside the Coulomb gap corresponding to equa- 
tion (|4]i. A super-linear conductivity power law was pre- 
viously observed in Si:As and Si:P by Castner and colla- 
borators II20II2 1 1 using resonator techniques at certain fre- 
quency within the range of the present work. Our results 
are in accord with the measurements on similar samples at 
higher frequencies (30 GHz to 3 THz) using optical tech- 
niques Il22li231 . In contrast, a sub-linear frequency depen- 
dence in the zero-phonon regime has been reported by Lee 
and Stutzmann [[161 based on experiments on Si:B in the 
microwave range and by Helgren et al. ||9l using quasi- 
optical experiments. 

3.2 Temperature dependence of the conductivity 

Leaving the zero-phonon regime by raising temperature, 
a gradual increase of the conductivity <ti is observed for 
all investigated Si:P-samples. First, the temperature depen- 
dence is approximately linear in agreement with the pre- 
diction of Austin and Mott |rTl,'24'|. With increasing T, it 
gradually becomes stronger until the charge carrier activa- 
tion into the conduction band dominates over the hopping 
transport. The onset of the thermal effects depends on the 
phosphorus concentration: the higher the doping, the lower 
the temperature at which the temperature dependence sets 
in. But the way it happens is the same for all the samples 
measured. Taking the example of n/ric = 0.83, we have 
plotted the conductivity values in Fig|5] at fixed frequen- 
cies of 150 MHz and 1 GHz. Above T = 9 K the curves 
merge into the thermally activated curve. 

In fig. |6] the frequency dependent conductivity is plot- 
ted for different temperatures, for the example of the crys- 
tal with the high concentration n/uc = 0.9; the conduc- 
tivity power law gradually decreases with rising tempera- 
ture. The transition to a sub-linear power law upon raising 
T is in accord with previous observations Il21ll24ll25l ; to 
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Figure 6 Temperature dependence of the conductivity 
spectra typical for all the investigated Si:P-samples. 
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Figure 7 Doping dependence of the low-temperature di- 
electric constant of Si:P. The solid line is the fit with the 
power law function ^ for the electronic dielectric func- 
tion El — £si- 



our knowledge, the theoretical description of the gradual 
decrease of a is lacking. 

3.3 Dielectric function The dielectric function si is 
independent of frequency in our frequency range, taking 
the measurement uncertainty into account. A fit with the 
function ^ results in an exponent ( =0.8, as shown in 
Figlll 

4 Conclusion Si:P-samples with relative doping n/ric 
from 56% to 90% show temperature-independent dynamic 
response at the temperature of 1.1 K. Over a wide fre- 
quency range (100 MHz to 5 GHz) the conductivity power 
law is slightly super-linear, with a faint decrease of the 
power a upon doping. At higher doping, above 90 %, the 
conductivity power law drops abruptly. Interpretation and 
further studies at lower base temperature using ^^He are in 



progress. At elevated temperatures the conductivity values 
are increasing and the power a is falling gradually with 
rising T. The electronic contribution to the dielectric func- 
tion is constant and increases with phosphorus concentra- 
tion as expected when the MIT is approached. 
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